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ABSTRACT 

We have generated accurate V and / template light curves using a combination of Fourier decom- 
position and principal component analysis for a large sample of Cepheid light curves. Unlike previous 
studies, we include short period Cepheids and stars pulsating in the first overtone mode in our anal- 
ysis. Extensive Monte Carlo simulations show that our templates can be used to precisely measure 
Cepheid magnitudes and periods, even in cases where there are few observational epochs. These tem- 
plates are ideal for characterizing serendipitously discovered Cepheids and can be used in conjunction 
with surveys such as Pan-Starrs and LSST where the observational sampling may not be optimized 
for Cepheids. 

Subject headings: Cepheids — distance scale — Cepheids — stars 



1. INTRODUCTION 

The Cepheid Period-Lmxiinosity (PL) relation is a fun- 
damental rung in the astronomical distance ladder. With 
the Hubble Space Telescope (HST) routinely resolving 
stellar populations in nearby galaxies, Cepheid distances 
can potentially be calculated for a large number of new 
systems. In this paper, we show that high quality tem- 
plate light curves can accurately fit pulsation periods and 
mean luminosities to sparsely sampled Cepheid observa- 
tions. It is our hope that these templates can be used to 
characterize Cepheids in archival as well as new observa- 
tions, even when the time sampling of observations has 
not been optimized for Cepheid studies. 

Principal component analysis (PCA) has found ap- 
plication in a wide range of astronomical studies (e.g., 
Il^abeil 119731 IConselicd ^OOl ) . The motivation behind 
PCA is that it can greatly reduce the number of vari- 
ables required to describe a data set. By performing an 
eigenvalue decomposition, PCA can be used to generate 
a small number of eigenvectors that describe the major- 
ity of the variance in a data set. In the case of Cepheid 
variables, one could use PCA analysis to construct light 
curves with the smallest possible number of free param- 
eters. Ideally, all Cepheids with identical periods would 
have a single unique light curve shape. In this case only 
a period, an average magnitude and a phase would be 
needed as free parameters to fit a light curve. 

The classical technique for determining a stellar 

pulsation period is to calculate a "stri ng length" 

(jLafler fc Kiil^^ [19651 iBiirke et all [l97(l . A hght 
curve is folded along a trial period and a string length 
is computed by summing distances between points con- 
secutive in phase. The trial period with the shortest 
computed string length is taken as the true period. This 
method is accurate for well sampled light curves with 
small photometric errors, but has been supplanted by 
recent techniques. 

More robust periods can be obtained by fitting tem- 
plate light curves. Fourier analysis of Cepheids was intro- 
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duced in lSchaltenbrand fc TammannI ()1971[ ). and lStetsonl 
(|1996f ) generated templates ba sed on Fourie r decompo- 
sition of MW and LMC stars. iTanvir et all ([2005,) used 
principal component analysis (PCA) to construct a large 
set of well sampled Cepheid obser vations. In this pa- 
per, we extend the techniques of T anvir et al.l (|2005( ) to 
include short period Cepheids as well as first overtone 
Cepheids. 

Short period Cepheids (with periods less than ^ 10 
days) have been excluded from many studies for a vari- 
ety of reasons. First, LMC Cepheids may show a dis- 
continuity in the PL-relation at 10 days, casting un- 
certainty on their utilit y as reliable standard candles 

geow fc Kanbuil[2008f ). Second, there is the possibil- 
ity of confusion between fundamental and first-overtone 
mode Cepheids that have similar periods but different 
PL-relations. Third, if fainter stars like short period 
Cepheids are used to estimate distances, incompleteness 
bias can skew the derived distances to smaller values 
(jSandage fc CarlsonI [l98l iFreedman et al.l [200 1[) . Fi- 
nally, the shorter period Cepheids show more variation 
in light curve shape, making template construction more 
daunting. Keeping these potential problems in mind, we 
boldly go forward and derive short period templates re- 
gardless. 

The outline of the paper is as follows: In Section [21 
we describe generating template light curves using PCA 
from a large literature sample of Cepheid stars. In Sec- 
tion [3] we perform Monte Carlo simulations to determine 
how precisely we can recover Cepheid parameters using 
our templates. In Section [H we discuss the process of 
converting fit parameters into a distance measure. 

2. PRINCIPAL COMPONENT ANALYSIS 

2.1. Training Set Selection 

Following iTanvir et al.l (|2005| ). we perform PCA on 
V and / band light curves simultaneously. We there- 
fore gathered light curves for stars that have both V 
and / measurements. The majority of our template 
stars come from the Optical Gravitation al Lensing Ex- 
perime nt (OGLE) databases for the LM C ([Udalski et al.l 
Il999a) and SMC ( Udalski et al.lll999bD. We g athered ad- 
ditional LMC observations from ISebo et all (120021) and 
additional LMC and SMC light curves from lMoffett et"aLl 
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Fig. 1. — Period distribution of the stars that are used to gener- 
ate our templates. The Milky Way, LMC, and SMC fundamental 
mode stars along with the combined LMC and SMC first overtone 
pulsators are plotted. 



(|1998f ). Our Galactic Cepheid sample was c ompil ed from 
light curves in the VizieR datab ase ([Oc hscnbci n etldl 
20001). including data prese nted in lBerdnikov, (|1997l ) and 
Berdnikov &: Turnerl ()2001[ ). We also included Galactic 
Cepheids from the McMaster Cepheid Photometry and 
Radial Velocity Dat a Archive wh ich includes light curves 
from many sources ( Gicrcn TgsT; Mo ffett fc BarnesHlQp: 
Coulson fc Caldwell .1985; Berdnikov fc Turner] Il995: 
Hendenlll996l : lBarnes et al.lll997D . We would have liked 
to include light curves from the MACHO survey. How- 
ever the MACHO survey uses non-traditional "blue" and 
"red" filters that are non-trivial to convert to V and /. 

To be included in our analysis we required each light 
curve have at least 15 epochs of observations in both 
V a nd I, with the exception of well spaced light curves 
from lBerdnikov fc Turner (2001), where we only demand 
6 epochs in both V and /. The initial sample contained 
3173 light curves, with 248 Milky Way LCs, 1378 LMC 
LCs, and 1514 SMC LCs. Unlike the majority of previous 
Cepheid studies, we did not exclude stars with periods 
shorter than 10 days from the analysis. We fit Fourier 
components to every light curve and rejected 305 out- 
liers (i.e., those with unusual light curve shapes or poor 
fits). The final sample included 150 MW Cepheids 677 
LMC Cepheids, 689 SMC Cepheids, and 1171 overtone 
pulsators from the LMC (549 stars), SMC (575 stars), 
and MW (47). The period distributions of these stars 
are plotted in Figure [TJ 

The OGLE database does not distinguish between fun- 
damental mode Cepheids and Cepheids pulsating in the 
first overtone mode. Figure [3] shows our simple cuts 
in period-magnitude space used to provisionally classify 
OGLE stars as fundamental or overtone pulsators. For 
the rest of the stars, we use the source catalog designation 
of overtone or fundamental Cepheid. Color-magnitude 
diagrams of the OGLE are plotted in Figure [2] 

We initially decompose each light curve in our sample 
into Fourier components by fitting the following equa- 
tions to the V and / light curves simultaneously. 
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Fig. 2. — Color-magnitude plots for the OGLE stars, uncorrected 
for reddening. Points have been color-coded by pulsation period. 
The upper panels show LMC stars while the lower panels show 
SMC stars. Fundamental mode pulsators are shown on the left 
and first-overtones are on the right. 

mv{t) = mv.o + E ^3 i'^^jt/P) + Pj cos (27rjt/P52) 

where P is the period, t is the Julian day of the ob- 
servation, m/.o a-iid mvfi are the average magnitudes in 
each band, and the a and (3 terms are the Fourier ampli- 
tudes. We constrain f3j=i to zero to impose a common 
phase for all the fits. These Fourier fits generate the 32 a 
and 13 values for each star that are used in the principal 
component analysis. While the shapes of the light curves 
are fitted simultaneously, the average magnitudes in each 
band are completely independent. This ensures our light 
curve templates are not dependent on the various dust 
corrections (or lack thereof) that have been applied in 
our different source catalogs. 

Ngcow ct al. (2003) discuss how fitting Fourier com- 
ponents to sparsely sampled data can result in poor fits 
(see their Figure 2b for an example of how Fourier decom- 
position can fail for sparsely sampled light curves). To 
keep our fits constrained to reasonable shapes, we create 
a smooth light curve by linearly interpolating each ob- 
served Cepheid to include 50 points evenly distributed 
across the full phase of the light curve. We then fit the 
Fourier components of this smoothed light curve, and use 
the results as the initial guess for the fitting of the ob- 
served data points. In the final fit, each Fourier compo- 
nent is allowed to change by a maximum of 20% from the 
smoothed fit parameters, thereby preventing divergences 
compared to the smooth fit. Figure S] shows typical re- 
sults for out fitted Fourier components. Figures [5] and [6] 
show the distribution of the 16 best fit a and [3 values for 
all of the fundamental and overtone /-band light curves. 

2.2. PCA Template Construction 

To test for differences in Cepheid populations, we con- 
structed templates based on various subsets of the data. 
In particular, we made templates that include all our 
fundamental mode Cepheids and subsets which included 
only LMC stars, only SMC stars, only short period stars, 
only long period stars, and only overtones. For each sub- 
sample, we first subtracted an average light curve from 
each measured light curve. We then performed a PCA on 
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Fig. 3. — First pass classification between fundamental mode and first overtone from the LMC and SMC OGLE data. Stars between the 
solid blue and dash-dot green line are provisionally labeled as fundamental mode and stars between the dash-dot green and dashed red are 
labeled as first-overtone. The lines have a slope of -3 and y-intercepts of 17.2,16.65,15.9 for the LMC and 17.8,17.15, and 16.5 for the SMC. 
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Fig. 4. — Example of an OGLE 2.8 day period Cepheid with the 
solid curves showing the best fit Fourier decomposition. The upper 
curve shows the /-band while the lower curve shows V. The sparse 
sampling in the V-band around a phase of 1.0 and 1.2 allows the 
fit slightly too much freedom (The curve rises slightly in those two 
regions). By performing the PC A analysis over many hundreds of 
light-curves, such discrepancies should be averaged out. 



TABLE 1 

Percentage of the total variance contained in 
each pca component for each of our 

constructed TEMPLATES. 



Model 


PCAl 


PCA2 


PCA3 


PCA4 


Total 




% 


% 


% 


% 


% 


short period 


70.1 


18.7 


4.8 


1.4 


95.1 


long period 


64.1 


16.6 


6.8 


2.9 


90.4 


LMC 


66.7 


15.6 


9.8 


1.8 


94.0 


SMC 


70.1 


17.2 


6.0 


1.5 


94.8 


overtones 


82.8 


6.0 


3.0 


1.0 


92.9 



the residuals to develop eigenvectors that can be used to 
correct for deviations from the average light curve. In all 
cases, over 90% of the variance could be described with 
only four principal component eigenvectors (Table [T]). An 
example of an average light curve and PCA eigenvectors 
are plotted in Figure [71 

We plot the PCA vector strengths as a function of pe- 
riod in Figurc[5] Several of the PCA vectors show strong 
trends with period. This is fortunate (and not too sur- 



prising), as it means while fitting for a star's period, we 
can simultaneously make an educated guess as to what 
the corresponding shape parameters should be. 

Before performing PCA, we needed 32 Fourier compo- 
nents, as well as two magnitudes, a period and a phase 
to accurately fit a Cepheid light curve. After performing 
PCA >90% of the variation in the light curve's shape can 
be described with just four eigenvector amplitudes. As 
a final step we note that the eigenvector amplitudes are 
strong functions of period, allowing us to make a quality 
template fit with only four free parameters (two average 
magnitudes, a phase, and a period). Figure [8] does show 
one possible pitfall as the first PCA vector for the short 
period stars shows a great deal of scatter and little trend 
with period. We therefore caution that it is possible this 
eigenvector amplitude should be left as a free parameter 
if possible when fitting a light curve. 

3. TEMPLATE FITTING ACCURACY 

iTanvir et al.l (|2005D have already demonstrated that 
the template fitting technique is superior to other com- 
mon period estimation techniques for cases where the 
photometry is noisy. They show magnitudes and peri- 
ods determined through template fitting can reduce the 
scatter in distance estimates by 30% compared to simple 
string length methods. 

We now endeavor to use Monte Carlo simulations to 
quantify how well our templates can recover magnitudes 
and periods from photometry of Cepheid stars. We 
have developed a fitting routine that uses Levenberg- 
Marquardt least-squares minimization to find the best 
fitting Cepheid light curve, given a set of V and / pho- 
tometry. The best fitting template is found by varying 
period, phase, < I >, and < V >. We have included 
an option to vary the amplitude of the PCA eigenvec- 
tors. Fitting sparsely sampled light curves, there is a risk 
of aliasing or converging on local minima. To avoid 
such problems, we use a series of initial guess periods and 
phases to ensure we find the global minimum. 

Our Monte Carlo varies four parameters to judge their 
impact on our fitting routine's robustness: First, we com- 
pare 5 well sampled Cepheids from the OGLE database 
with different periods; Second, we vary the total number 
of observations in each band; Next, we look at possible 
effects of template mis-match (fitting LMC stars with a 
template derived from SMC stars, fitting overtone stars 
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Fig. 5. — The distribution of Fourier decomposition parameters for the best fits to all fundamental mode Cepheids in the /-band. For 
each population of stars, we fit low order polynomials to these Fourier component versus period plots and reject stars from further analysis 
that are more than Scr outliers. The number of stars in our sample is large enough that we have plotted contours in the densely populated 
regions and individual points in the sparsely sampled regions. 



TABLE 2 

The different epoch samplings used in our 
Monte Carlo error analysis. 



Model 


/ observations 


V observations 


total 


1 


20 


15 


35 


2 


10 


5 


15 


3 


5 


3 


8 


4 


4 


2 


6 



with fundamental mode templates and vice- versa); Fi- 
nally, we vary the photometric precision of the light curve 
points. 

For each realization the star was randomly sampled in 
both bands. Therefore there are Nj + Ny unique epochs 
of observations and we do not explicitly model observing 
strategies that observe in multiple filters simultaneously. 

We used four fundamental mode Cepheids (periods of 
2.5, 4.2, 16.0, and 20.7 days), and one overtone (period 
2.5 days) from the OGLE LMC data set for the Monte 
Carlo tests. All the photometry has initial errors of or- 
der 0.013 mags in both bands. Although there could be 
some objection to using stars which are included in the 
PCA analysis, our sample sizes are large enough that the 



addition or subtraction of a few stars should make little 
difference to our final templates. 

For the error analysis, we fit for only period, my, mi, 
and phase. We rely on the polynomial fits in Figure [5] to 
give reasonable amplitudes for the PCA eigenvectors. We 
also folded the original light curves so that the photom- 
etry covers at most 3 periods. This restriction is needed 
to keep the fitting procedure from falling into local min- 
ima caused by aliasing. When fitting sparsely sampled 
data over a long baseline a more brute-force exploration 
of period parameter space would be required than our 
current fitting procedure. 

Figure [11] shows an example of one of our fitting simu- 
lations. As expected, the errors are largest for the over- 
tone Cepheid when fit with the wrong template. We also 
show the calculated uncertainties reported by our fitting 
routine in Figure [TTl In general, the reported uncertain- 
ties are a good match to the actual errors resulting from 
the fits. The one exception is that the overtone uncer- 
tainties are under-estimated as a result of making the 
assumption that the reduced should be unity, which 
is clearly incorrect in the case of fitting an overtone with 
a fundamental mode template. When we use an overtone 
template, there is a clear improvement in the values, 
indicating that it is a better fit. 
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Fig. 6. — Best fit Fourier components for the overtone Cepheids in the /-band. The distributions for the high order coefficients 
are concentrated near zero, implying these light curves are simply sinusoidal, an unsurprising result for stars commonly referred to as 
"sinusoidal" or sCepheids. 



We have also explored possible template mismatches 
due to metallicity by repeating the Monte Carlo experi- 
ment using an SMC template ([Fe/H]- -0.7) to fit LMC 
stars ([Fe/H]'^ —0.3). The resulting periods and magni- 
tudes are practically identical, suggesting that the tem- 
plates can be used across different metallicity popula- 
tions. 

A proper assessment of the errors of individual mea- 
surements is key for determining a proper P-L relation. 
We have therefore compared the uncertainties reported 
by our fitting procedure to the true offsets seen in our 
Monte Carlo re-sampling. Overall, the returned uncer- 
tainties are comparable to the actual errors derived from 
the Monte Carlo analysis. When the photometric er- 
rors were low (0.01 mag), the returned uncertainties were 
slightly too small (by ~15%). On the other hand, when 
the errors were high and the sampling was sparse, the 
returned uncertainties were slightly larger than the true 
Monte Carlo calculated errors. Typically, around 10% 
of the fits would fail catastrophically, usually caused by 
aliasing or very sparse sampling. These failures could 
readily be seen as poor values or by visual inspection. 
This seems to imply that we can use the returned uncer- 
tainties, but as we will see later, observing strategy also 
affects the results and it is probably best to simply run 
a Monte Carlo for the photometric error and uncertain- 



ties for the specific sampling frequency used by a given 
observing program. 

We ran several simulations to test how sensitive our fits 
are to the phase coverage of the observations. Generally, 
if the observations do not span more than half of the 
full phase, there is a large likelihood (10%-40%) that the 
fitted period will catastrophically fail (defined as a final 
fitted distance error of >10%). 

To illustrate how well our fitting procedure works, we 
compared our template fits to fits using a simple asym- 
metric saw-tooth function. The results are plotted in 
Figure [131 When the light curves are well sampled, the 
templates and saw-tooth converge to practically identi- 
cal values. In the sparsely sampled case (only 7 observa- 
tional epochs) , the templates return accurate fits in 4 out 
of 5 cases while the saw-tooth function fits fail in every 
case. 

In summary, our tests indicate that our PCA template 
technique can fit periods with an precision of ±0.1 — 0.3 
days with only 6 epochs and photometric precision of 
0.01 mag. If the number of epochs increases to 15 days, 
our uncertainty drops to ±0.03 — 0.2 days. 

4. CONVERTING FITS TO DISTANCES 

Having established that our templates can accurately 
fit Cepheid periods and average magnitudes, we now 



6 



Yoachim et al. 




Phase Phase 

Fig. 7. — The average light curve (a) and eigenvectors from the fundamental Cepheid principal component analysis, /-band is shown 
in thick red while V-band is shown in blue. The first PCA vector (shown in b) primarily controls the magnitude of pulsation, while the 
other PCA vectors (c-e) control the exact shape of the sawtooth rise of the light curve. Template light curves are constructed as linear 
combinations of these vectors, an example of which is shown in (f). Typical coefficients for these vectors are plotted in Figure [8] 
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Fig. 8. — The template PCA coefficients (all multiplied be a factor of 10^) for all the fundamental mode Cepheids. The SMC, LMC, 
and MW stars are shown by red, blue, and black points respectively. Our polynomial fits are plotted as solid curves. We have separated 
long period, short period, and overtone Cepheids into the left, middle, and right panels respectively. 
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Fig. 9. — Examples of four OGLE fundamental mode Cepheids fit with our derived templates. Both the V and I light curves are fit 
simultaneously. Dotted lines show the best-fitting curve if we fit with just the average light-curve (PCA eigenvectors fixed at zero), dashed 
curves show the best fit if we set the PCA eigenvectors to the best polynomial fit values for the given period, solid curves show the fit if 
we leave the first two PCA vectors as free parameters. The numbers in parentheses show the total number of free parameters in each fit. 
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Fig. 10. — Same as Figure |9] only now for Cepheids identified as first-overtones and fit with our overtone template. Again, it does not 
seem necessary to leave the PCA amplitudes as free parameters to converge on a quality fit. 



point out some of the potential pitfalls in using fitted 
parameters to derive accurate distances. In theory, dis- 
tances from single Cepheids can be averaged together to 
find a precise distance to a galaxy. Besides the usual 
systematic errors associated with photometry, there are 
additional caveats that apply when observing a sample 
of Cepheids: (1) Cepheids used for a distance calcula- 
tion must be above the completeness limit of the obser- 
vations, otherwise, faint short-period stars will be under- 
sampled and the final calculated distance will be biased 
(|Sandage fc CarlsonI [19881 iFreedman et al.l[2?)0l : (2) If 
the Cepheids do not sample the full range of the instabil- 
ity strip, the y can be offset from a standard PL-relation. 
iMager et al.l (IIOO8) discuss how the scatter in the PL- 
relation in the outskirts of NGC 4258 (|Macri et al.ll2006D 
is greatly reduced because the Cepheids populate a lim- 
ited region of the instability strip; (3) Finally, there is 
always the risk that a Cepheid may not be de-blended 
from a nearby optical/physical companion. Blending 
is expected to bi as Cepheid distance measurements to 
smaller val ues (fStanek fc Udalskilll999l : iMochejska et al.l 

l2000l[200l . 

In a companion paper (|McCommas et all [20091 ) . we 
have successfully applied our templates to HST data. 
Figure [TJj shows some examples of how well the tem- 
plates can fit noisy and sub-optimally sampled data. The 
templates only start to fail for the star with the longest 



period, where only half of the full phase is observed. 

4.1. Metallicity Effects on Light Curve Shape 

There is some question as to whether the shape of 
the light curve can be used to determine t he metallicity 
of a Cepheid (iPaczvnski fc Pindorl 120001 : iKanbur et al.l 
l2002f ). Looking at the long period Cepheids (Figure [8]), 
the first and second principal components have system- 
atically smaller values for the SMC stars compared to the 
LMC and MW, suggesting the shapes might be intrinsi- 
cally different. We use a 2D Kolmogorov-Smirnov test 
to compare the distributions in the PCA-Period distri- 
butions (the MW, LMC, and SMC long period variables 
all have very similar period distributions). The LMC 
and MW are consistent (Probability > 0.12) with being 
drawn from the same populations for all of the PCA vec- 
tors (i.e., the shape of the light curve at a given period 
does not show significant differences). The SMC shows 
a significant difference between the MW (P=0.008) and 
LMC (P=0.04) in the distribution of the first PCA vec- 
tor, but not the higher order PCA vectors. While the 
SMC light curve shape is different on average, there is 
little hope of assigning a metallicity to an individual star 
based on its light curve shape, since the intrinsic scat- 
ter within the SMC is larger than the differences from 
the LMC or MW. Even with the large number of stars 
in our sample, it is difficult to differentiate between the 
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Fig. 11. — Results of our Monte Carlo simulation where we recover Cepheid properties using our templates and least-squares fitting 
procedure. We fit 5 different stars, along with the 5 different observing strategies listed in Table [2] The four panels on the left show the 
errors in fitting magnitudes and periods and the corresponding error in distance modulus. On the right, we show the uncertainties returned 
by our least-squares fitting routine. With photometric errors of 0.01 mags and only 6 total epochs of observations, the uncertainties of the 
fitted parameters result in only a 0.1 mag uncertainty in Afi'^^^. 



low metallicity SMC Cepheids and the higher metallic- 
ity LMC and MW Cepheids. It is possible, in theory, to 
observe enough long period stars that one could distin- 
guish if a population was high or low metallicity based 
on the PCA distribution. However, it would require a 
prohibitively large sample size, as well as sufficient phase 
coverage that the shape parameters could be measured 
robustly. 

Unlike the long period variables, the short period 
Cepheids in the three systems have very different period 
distributions. However, the PCA distributions overlap 



to a sufficiently large degree (Figure [8]) that light curve 
shape cannot readily be used to measure the metallicity 
of individual stars. 

While we find metallicity does not alter light curve 
shapes significantly, there is evidence that metallicity dif- 
ferences can alter the PL-relation . Several stu dies claim 
to find a nietallic it y dependence (lKennicutt~e t al. 1998; 
'Sakai et al .j Eool iMacri et a l.' '2006^, 'Saha et al. 2006; 
Romaniel lo et al.l 120081 : ISandage fc Tammann 2008,), 
while ot hers find that the PL-relation is constan t across 
galaxies (jUdalski et al.ll200lHGreren et al.ll2005[ ). 
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Fig. 12. — Same as Figure [TT1 only now the photometric errors have been increased to 0.2 mag. Again, the actual errors (left side) and 
very similar to the returned uncertainties (right side). 



5. CONCLUSIONS 

We have e xtended the princip al component analysis 
techniques of iTanvir et al.l (|2005| ) to generate template 
light curves of Cepheid variables and first-overtone vari- 
ables. We have used a Monte Carlo simulation to demon- 
strate how robustly our templates can be used to fit 
Cepheid periods and magnitudes. Unlike previous stud- 
ies, we do not limit ourselves to stars with periods longer 
than 10 days. Finally, we demonstrate the effectiveness 
of our templates on HST data and show that our tech- 
niques can be used to fit accurate periods and luminosi- 
ties even when the observations have not been optimally 
spaced for observing variable stars. Our templates open 



up a new regime of distance measurement possibilities 
by enabling accurate fits for long period, short period, 
and overtone Cepheids from noisy and sparsely sampled 
observations. 
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Fig. 13. — An example of how effective our templates can be at fitting sparsely sampled data. In the left hand panels are the /-band light 
curves from Figure |9] fitted with our LMC templates (blue) and a simple asymmetric saw-tooth function (red). The right hand panel shows 
the results when 7 random points from the light curve are sampled and the the fitting is repeated. In four of the 5 cases, the templates 
are able to fit accurate average magnitudes and periods. The larger number of free parameters in the saw-tooth function prevents it from 
properly converging in the sparsely sampled cases. 
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Fig . 14. — Example of V and / light curves of Cepheids in M81 observed with HST and fit with our templates. See lMcCommas et al.l 
II2009I ) for a larger sample of Cepheids and a full distance calculation to M81. The fitted period in the lower right is poorly constrained 
due to the sparse phase coverage. 
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